Zeeman smearing of the Coulomb Blockade 
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Charge fluctuations of a large quantum dot coupled to a two-dimensional lead via a single-mode 
good Quantum Point Contact (QPC) and capacitively coupled to a back-gate, are investigated in 
presence of a parallel magnetic field. The Zeeman term induces an asymmetry between transmission 
probabilities for the spin-up and spin-down channels at the QPC, producing noticeable effects on the 
quantization of the grain charge already at low magnetic fields. Performing a quantitative analysis, 
I show that the capacitance between the gate and the lead exhibits — instead of a logarithmic 
singularity — a reduced peak as a function of gate voltage. Experimental applicability is discussed. 
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Recent research on mesoscopic quantum dots has led 
to a revival of Kondo physics. There is an extensive liter- 
ature on the (one-channel) Kondo behavior of small dots 
with an odd number of electrons and a finite level spac- 
ing e, attached to two leads. In such an arrangement, the 
dot acts as an Anderson impurity jlj. In addition, it was 
shown by Matveev g that the Hamiltonian of a large dot 
(e — * 0) in the Coulomb blockade regime, coupled to a 
two-dimensional electron gas (2DEG) via a good single- 
mode Quantum Point Contact (QPC) and capacitively 
coupled to a back-gate, is mathematically equivalent to 
the two-channel Kondo model || in the Emery-Kivelson 
limit H (Fig. 1). The two channels of the Kondo prob- 
lem correspond to the two spin channels for tunneling 
through the QPC. For recent reviews, see Refs. ||||. 

The non-Fermi-liquid nature of the ground state of the 
two-channel Kondo model is here reflected by the non- 
analyticity of the capacitance measured between the gate 
and the reservoir near the points where the dot charge Q 
is half-integer § (C = dQ/dV G ) 

C(N) =C - bC gd \X\ 2 cos( 27 rA0 In ^ ■ (1) 

N is a parameter proportional to the gate voltage eN — 
VaCgd, C g d denotes the gate-dot capacitance, C Q is the 
total dot capacitance and |A| 2 <§C 1 is the small reflection 
probability at the QPC. Here, b > is proportional to 
the Euler constant 7 = e c , with C sa 0.5772. The sec- 
ond term describes the cross-over from the linear charge- 
voltage dependence to the "Coulomb staircase" behav- 
ior (inset in Fig. 3). This effect has been observed by 
Berman et at, in AlGaAs/GaAs heterostructures |@,|[. 

Below, I investigate how the capacitance C(N) evolves 
if a magnetic field B is applied parallel to the 2DEG 
||. I will show that the logarithmic divergence in the 
differential capacitance 6C(N) = C(N) — C with N — > 
1/2 gets already cutoff by a small magnetic field (Fig. 3) 

^-(N) =6|A| 2 (l- ( 5 2 )cos(2 7 rA^) x (2) 

Cgd 

lnMax[|A| 2 <5 2 sin 2 (7riV); |A| 2 cos 2 (ttA0] , 



< 5 -C 1, being proportional to the Zeeman energy, 
A = gfisB. This results from a Zeeman-like asymmetry 
between reflection probabilities for the spin-up and spin- 
down channels at the QPC. For a high field (5 as 1), only 
the spin-up channel will be transmitted producing a one- 
channel QPC model. The capacitance (or charge Q) only 
exhibits periodic oscillations as a function of N |2j , 

C(N) = C - bC gd \X\ cos(2irN). (3) 

The Coulomb staircase behavior gets completely smeared 
out (inset in Fig. 3). An experiment using AlGaAs/GaAs 
heterostructures is presumably appropriate to probe this 
effect in the capacitance jlo| . Indeed, Zeeman splitting at 
the QPC in an in-plane magnetic field has been confirmed 
by Thomas et al via conductance measurements [fllf , and 
for few conducting modes at the QPC the Lande factor 
is enhanced g ~ 1 |Q (the bulk value is \g\ — 0.4). 

Let me emphasize that below a quantitative analysis 
of the smearing out of the logarithmic peak for the ca- 
pacitance 6C(N) is performed (which has not been pre- 
viously done in Ref . j^j ) . The crossover from the two- to 
the one-channel QPC model is carefully investigated. 

First, it is useful to compute reflection probabilities 
at the QPC in presence of a magnetic field parallel to 
the 2DEG, and to discuss the necessary magnetic-field 
dependent adjustment of the QPC potential V a (B). In 
the close vicinity of the contact, the (smooth) confining 
potential will be approximated as a harmonic one H . 
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FIG. 1. The experimental setup: A large dot coupled to a 
2DEG via a single mode QPC and capacitively coupled to a 
back-gate, in a parallel magnetic field B (C a = C s d). 
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An electron with spin-projection a =t,4 (or a = ±1) 
along the magnetic field axis, is then subjected to the 
total potential V a (B) — muj x 2 x 2 /2 + mL0 y 2 y 2 /2 — aA/2. 
Here, u x and ui y are the curvatures of the potential, and 
m the mass of an electron. The transverse part of the 
Hamiltonian produces "n transverse modes" . Then, the 
one-dimensional (ID) wave function ^(x) for motion 
along x is determined by the effective potential 

V£(x) = V a (B) -a^ + hw y (n + ~) - \mu x 2 x\ (4) 

the explicit form of which is of no interest here. The 
height of the barrier potential at the QPC for a spin- 
down electron is enhanced by A. From now on, 1 ad- 
just V a (B) =V a (B = Q) + A/2, because I want field- 
independent reflection probabilities for spin-up channels. 

Near the saddle-point (x ~ 0), the threshold energies 
of the mode n are then spin-splitted, as follows 

El=E n = V a (B = 0) + hu v (n + l), (5) 
E 7 [ = E n + A. 

Classically, modes with threshold energy below the Fermi 
energy Ep are perfectly open and the others remain 
closed. But, quantum mechanically transmission and re- 
flection at the saddle arc neither completely open nor 
completely closed [O]. Here, I fix the voltage V a (B = 0) 
such that Eq <C Ep < E\. In Fig. 2, this corresponds to 
adjust £ such that 0.5 <C £ < 1.5 (lu v /uj x = 1). Then, the 
spin-up channel of the transverse mode n=0 remains at 
almost perfect transmission whatever the applied mag- 
netic field. Moreover, I can disregard modes n > 1 which 
are almost perfectly reflected, because I am interested 
only in transport through the constriction. 

The reflection/transmission amplitude of a ID particle 
passing through an inverted parabolic barrier has been 
studied in detail by Connor H|. Taking n = 0, small 
reflection probability for the spin-up channel reads 

Kl = \\\ 2 = - 1 —— - ps exp(-2^ T ), (6) 

1 + cxp(2-7t£q ) 

where £q = [Ep — Eq]/Tiuj x ^S> 0. Below, I will focus 
on low magnetic field effects. Using Eq. (|J), similarly I 
obtain (5 = A/2Huj x < 0.4) 

K l ps cxp(-27r5 i ) ps nl(l + 471-6). (7) 

Both channels are transmitted but 1 3> TZq > T^J (See 
Fig. 2). Applying an in-plane magnetic field, one gets 
what I call a two-channel anisotropic QPC model; This is 
defined, below. The limit of strong fields will be reached 
when the Zeeman energy approaches the curvature ener- 
gies of the potential: = E x and = 1 - ps 0. A 
single channel will subsist in the constriction. 
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FIG. 2. Exact transmission probability for the transverse 
mode n = and the spin-channel opposite to the magnetic 
field, as a function of £ = [Ef — V a (B = 0)]/hui x . The ratio 
ijj y /uj x , being 1. The different curves correspond to different 
values of the magnetic field i.e. 8 = A/2Huj x = 0, 0.1, 0.2, 0.4. 

As I am interested in the dynamics of the system at 
energies much smaller than Ef, I may linearize the spec- 
trum of the lD-fermions in state ^ a {x) = 4'™ =0 (a;). 
One can always write ^> a (x) = exp^kpx)^ R a (x) + 
exp(— ikpx)^! La{x); ^ La and ^R a respectively describe 
left- and right moving fermions, and Ep — kp 2 /2m. Fi- 
nite reflection in the channel a =t, J. then can be simply 
accounted for by adding a backscattering term f^||[l5| 

H bs = %VW*L(°)*fl«( )- (8) 

Vf denotes the Fermi velocity |l6) . One must also include 
the kinetic energy through the constriction, 

/+oo 
dx{^ Ra d x * Ra - W La d x ^lt a \- (9) 
-oo 

At almost perfect transmission, the electronic wave func- 
tion is shared between the reservoir and the dot. I can ne- 
glect finite size effects in a dot at the micron scale (e — > 0) 
§,m|l7). Note that (H km + H bs ) with 1 > TZ l > ftj de- 
scribes a two-channel anisotropic QPC model. Again, I 
ignore higher modes confined to the reservoir, and also 
neglect the Pauli contribution of the 2DEG. The charging 
process is described by the following usual term 

h c = e c {q-n)\ (10) 

with E c = e 2 / (2C g d) <§C Ep the (charging) energy that it 
costs to transfer a particle from the lead to the dot. The 
charge Q (of the dot) in H c is now normalized to e, 

r + oo 

Jo 

At low energies, I can proceed with this model by 
bosonizing the ID Fermi fields [^|,|(i| Jl5| , [l^ 1 . 




(12) 
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a is a short-distance cutoff, again a = ± for spin up and 
spin down, and p = ± for right and left movers. The 
spectrum of ID free electrons yields separation of spin 
and charge. Resulting Hamiltonians are plasmon-like 



Etkin — 



j=c,s 



/ + OG 
dx [{d x 
-co 



/) 2 + n/]. 



(13) 



d x <j)j with j = (c, s) measures fluctuations of charge/spin 
density in the constriction and IF, = d x 9j being its con- 
jugate momentum. In this representation, the charging 
Hamiltonian H c reads p,p||lq] 



H C = E C [^ 



(14) 



To minimize energy, the charge in the dot is pinned at 
the classical value Q c = (f> c (0)^/2/^ w N §,|@. 

Now, one has to examine the quantum corrections to 
the charge entering the constriction, in presence of a 
small (parallel) magnetic field. For energies smaller than 
E c and |A| <1, I can replace cos[\/27r0 c (O)] by the aver- 
aged value yJ^E c a/nv F cos(7riV) (similarly for the sinus). 
The prefactor comes from the average of the cosine term 
over the ground state of Hkm 20 |l5|]. Using Eqs. (p|),(Q) 
and ( p^ ) the backscattering term then reads 

|A| 



H bs » ^^-fv F E c acos(TrN) cos[V4^ s (0)] (15) 
Zira 

\X\S 



+ — Jjv F E c a sin(7riV) sin[V47r$ ^ (0)] . 

The previous anisotropy parameter 5 <C 1 is equal to n8. 
As in Ref. Q (page 417), I introduce the symmetric com- 
bination of the spin Bosc fields with respect to the QPC, 
(j) s {x) = [$ s (x) + $„(-x)]/y/2. The kinetic part is plf 



H kin [^ s 



v F 



+ OC 



dx (d x $ s y 



(16) 



In the absence of the magnetic field, i.e. 5 = 0, the 
Hamiltonian Hki n [&s] + H^s is equivalent to the two- 
channel Kondo model at the Emery-Kivelson line 



HeK = Etkin[^s] 



J± X [N] 



tt a 



cos[V47r$ s (0)]5 x , (17) 



with the Kondo exchange J± x = \\\^f^v F E c a cos(irN). 
S x describes the x-component of a fictitious quantum 
impurity. This equivalence occurs because S x commutes 
with H E k i). One can set S x = 1/2 (or S x = -1/2 &). 
H F k can be exactly refermionized and diagonalized []4[g| . 
The backscattering contribution to the ground state en- 
ergy then takes the form SE(N) = -~^-\n(E c /T(N)), 
where T(iV) = J± x 2 /irvpa is the Kondo resonance. The 
quantum correction SQ to the charge in the dot, becom- 
ing equal to SQ = -dSE/(E c dN) |J|. This results 
explicitly in SQ w -bill sin(27riV) \n{E c /T(N)) f§. 




FIG. 3. Normalized differential gate-lead capacitance as a 
function of N, for different values of the anisotropy parameter 
5 = nS = 0, 0.1, 0.2, 0.4, 1 (from top to bottom). In Fig. 2, I 
fix £ = 1.15 such that 72.J = 0.25, and 6 = 1. I also draw the 
total charge Q — Q c + SQ as a function of N for 5 = 0, 1. 

The charge fluctuations get "strongly" influenced by 
the spin ones; The logarithmic factor in SQ comes from 
spin fluctuations at the QPC For 5 = 0, these get 

pinned at the time h/T at which the screening of S x takes 
place and the Kondo coupling J± x diverges. For N=l/2, 
T = 0: There is a degeneracy between the dot states with 
Q=0 and Q=l [See Fig. 3; Q = Q c + SQ}. The capaci- 
tance displays a logarithmic divergence at this degener- 
acy point [Eq. (g)]. If JV < 1/2 one gets v^M ) ~ 
and V4^$ s (0) » 2tt if N > 1/2 [||. For 5 ^ the 
backscattering contributions to the ground state energy 
can now be calculated using the two-channel anisotropic 
Kondo model at the Emery-Kivelson line |24[] 



h ek = Hek 



sin[V47r$ ;j (0)]S* y 



(18) 



The extra magnetic exchange is (See Ref. M , page 394) 



(19) 



J_Ly[N] = \\\5^/~/v F E c asiii(TrN), 

S y being the y-component of the fictitious impurity. 

The equivalence between Hki n [$s]+Hbs and Hek can 
be understood, as follows. Starting with low-magnetic 
fields and at almost perfect transmission, the following 
commutators can be neglected whatever N G [0; 1] (for 
frequencies larger than u> m i n ; See discussion below) 



[H£ K ,§ x ]<xi&\\\§,<*:i 

[H§ K ,S y ](xi\X\S z <^l. 



(20) 



This implies that S x and S y both can be considered as 
good quantum numbers: The impurity spin in H^ K can 
oscillate between the two states S x = 1/2 and S y = —1/2. 
The use of the two-channel anisotropic Kondo model is 
completely justified, for all N G [0; 1]. It is noteworthy 
that this model can be refermionized and diagonalized; 
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See e.g. Ref. ||, page 401. Another energy scale re- 
lated to the magnetic field emerges: T = J± y /irvpa = 
-fE c lzl5 2 sin 2 (7riV)/7r. This is associated with the diver- 
gence of J± y . Using the two-channel anisotropic Kondo 

model, for B ^ I find (^sm{ln[iT-u]} = - taxT^/u]) 



5E = + — 9m{ Infer -oj] + hifeT - oj]}. 

Ju min 



(21) 



Let me comment on the low-frequency (spin) cutoff. 

Away from half-integer values of N, T tends to zero 
because sin(7rA) <C 1. Then, the coupling J± y < 1 is 
not really affected by the renormalization of the short- 
distance cutoff a. As in the absence of the magnetic 
field, all spin fluctuations get pinned at the time Ti/T. 
The equilibrium state at the QPC then corresponds to 
V4^$ s (0) w 7T [2tt] if N < 1/2 [if N > 1/2], and 
J±yS y < sin[V^7r$ s (0)] > w 0. Approaching the point 
N = 1/2, in contrast the Kondo resonance vanishes. But 
now, T becomes dominant such that J± y is flowing off 
to strong couplings at the finite time %/T. Spin fluc- 
tuations at the QPC then get pinned also for N ps 1/2 
and V4~7r<f> s (0) ps 37r/2. In the presence of the magnetic 
field, therefore the logarithmic integral on spin fluctua- 
tions in Eq. ( plj ) now acquires a natural low-energy cut- 
off, huijnin — Maxp; T], whatever the value of N £ [0; 1]. 
Performing the integration in Eq. (^l|), then I precisely 
obtain SQ ps -67^(1 - S 2 ) sin(27rA) hi(E c /huj min ). 

The charge (and spin) in the dot becomes a "continu- 
ous" function of N. 

This enables me to allege that the capacitance C = 
dQ/dVo displays, instead of a logarithmic singularity, a 
reduced peak as a function of N; See Eq. (||) and Fig. 3. 
For a small deviation n = (N — 1/2) from N = 1/2, one 
gets r(n) oc n 2 and T(<5) oc 5 2 , and then 



SC(x 



In 6 : n -C S 
Inn : ii<n. 



(22) 



Increasing the in-plane magnetic field, i.e. <5, then the 
capacitance peak reduces and becomes more and more 
broadened. For high fields 6 ps 1 and S ps 0.4, only 
spin-up electrons can be transmitted: Tq — > 0. One 
must disregard the spin-down channel at the QPC. For- 
mally, htu m in now exceeds -E c and the classical value of 
the grain charge must be rescaled as: Q c — 0|(O)/y / 7r = 
[</» c (0) + s (O)]/\/27r » N. Now, this produces £E = 
— |A|7^ C cos(27rA^)/7r 2 and then only a weak residual os- 
cillation in SC = -b\X\ cos(2tt7V) for all N £ [0; 1], as for 
"spinless" fermions j2j. The logarithmic part disappears. 

To conclude, investigating zero-temperature properties 
of a large dot coupled to a two-channel anisotropic QPC, 
I have shown that the steps of the Coulomb staircase 
function gets immediately smeared out by a finite mag- 
netic field parallel to the 2DEG. This rigorously proves 



that like the magnetic susceptibility of the Kondo impu- 
rity in the two-channel Kondo model jp4|, the logarith- 
mic singularity in the capacitance is destroyed by a small 
channel anisotropy. For high fields, adjusting correctly 
the voltage at the QPC one reproduces a single-channel 
QPC model; For |A| -C 1, the charge entering the dot 
exhibits only a weak quantum oscillation around its clas- 
sical value. Experimentally, the prominent smearing of 
the Coulomb staircase — predicted by increasing the in- 
plane magnetic field — could be possibly observed for suf- 
ficiently low temperatures T, probing (preferably) the ca- 
pacitance line shapes of the large dot in a single-terminal 
geometry (!(]] . To account for nonzero temperature, one 
must rescale Max[F; T] as ksT + Maxp; T]. Experimen- 
tally, it would be crucial to minimize the broadening of 
the capacitance peaks due to thermal effects. 

This work is supported by the Swiss National Science 
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